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Abstract 

In this paper we present a numerical scheme for the approxima- 
tion of singular shock solutions of the Keyfitz-Kranzer model system. 
Consistence in the sense of distributions is studied. As long as some 
numerical properties are verified when the space step tends to 0, we 
prove that the scheme provides a numerical solution that satisfies the 
equations in the sense of distributions with an approximation that 
tends to when /i — )• 0. We also show that this scheme adapts to 
degenerate systems. This is illustrated by two examples: the sys- 
tem presenting delta wave solutions originally studied by Korchinski 
and another system studied by Keyfitz-Kranzer that models elasticity. 
Consistence of the scheme in the sense of distributions is fully proved 
in the case of the Korchinski model. 

1. Introduction. Singular shocks have been discovered and investigated 
by different authors, see in particular [1,4,5,6]. They have been observed from 
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various viscosity techniques: Dafermos-Di Perna viscosity in [1] , usual viscos- 
ity in [4]. In the case of singular shocks viscosity solutions converge so weakly 
that their pointwise limits do not satisfy the classical Rankine-Hugoniot con- 
ditions. Besides this fact a unique entropic solution to the Riemann problem 
has been obtained in [1] for arbitrarily large data. In this paper we propose 
a numerical scheme based on a splitting technique that captures the singular 
shocks. We observe results exactly similar to those obtained in [1,4] with 
their respective viscosity techniques. Studies have shown the relevance of 
this scheme for other systems presenting irregular solutions. In our study of 
irregular shocks we consider two standard first order model systems of two 
equations whose solutions of the Riemann problem involve singular shocks 
and delta shocks. We also notice that this scheme provides neat results for 
the Keyfitz-Kranzer system of elasticity [2] for which the intrinsic difficulty 
is different from those in the two systems above. 

This paper focusses on the Keyfitz-Kranzer system 

ut + (m' - v)^ = 0, (1) 

vt + {\u^ - u)^ = 0, (2) 
which produces singular shocks [1,4], and the system 

Ut + {u\ = 0, (3) 

Vt + {uv)^ = 0, (4) 

originally considered by Korchinski [3] who discovered and investigated delta 
shocks in the solution of the Riemann problem. 

Let Uh, Vfi be the sequence of approximate solutions from the scheme. 
Under simple numerical properties to be rigorously proved, or to be admitted 
from numerical tests, we prove that the scheme is consistent in the sense of 
distributions in the following sense: V(0, V') e (C^(M x M+))^, 

Uh4)t+{{uhf-Vh)(f)x]dxdt^0, J[vhi't+{^{uhf-Uh)i^x]dxdt^O, (5) 

respectively J[uh(j)t + {{uhY)(j)x]dxdt 0, j[vhiJt + {uhVh)i>x]dxdt 0, 
when the space step /i — > 0. This means that the functions Uh,Vh tend 
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to satisfy the equations when — > 0. 



For system (1,2) we check numerically that the needed assumptions are 
satisfied for values of h as small as possible. We rigorously prove that, in 
the case of system (3,4), for any initial condition vP G L^(M) n L°^(M) and 
v° e -t'^(IR), these assumptions are satisfied. Therefore the scheme is consis- 
tent in the above sense. Of course, in the first case, from a rigorous point of 
view, one cannot be sure that these numerical assumptions always hold for 
every h when /i — )> 0. The proof in this paper shows that, for any given fam- 
ily of test functions with uniformly bounded support and uniformly bounded 
first and second derivatives, then a numerical solution satisfies the equations 
in the sense of distributions within a small deviation depending on h when- 
ever these assumptions remain valid. 

2. A numerical scheme. The singular shocks of the Keyfitz-Kranzer 
equations are unbounded which makes the elaboration of numerical schemes 
difficult: in the scheme below the numerical velocity u in system (1,2) can be 
unbounded when the space step h tends to which forces us to accept that 
the CFL coefficient r tends to when /i — > in order to preserve the CFL 
condition < 1. Therefore r = r/j depends on h and also on time so 

that r/i||wfe||Loo < 1. 

If Th tends to ( i.e. if tends to oo) slowly enough, then one can 

nevertheless obtain a convenient numerical scheme, although of an order less 
than one, on condition that for each iteration the assumptions are verified 
when /i — > 0. This ensures consistence of the scheme in the sense of distribu- 
tions, although the limit is not a distribution in general: it can be a singular 
shock in the case of the Keyfitz-Kranzer equations. Numerical results are 
given to prove that the set of assumptions is satisfied in representative sit- 
uations of singular shocks. In the case of the Keyfitz-Kranzer equations the 
scheme consists in a splitting of equations into the two subsystems 



vt + {vu)^ = 0, 



(6) 
(7) 



which is treated by transport with velocity it, and 



Ut = V, 



'XI 



(8) 



Vt — {vu 




XI 



(9) 
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which is treated by a centered discretization. In between, we introduce an 
average step in u, v which is needed in general to avoid oscillations due to 
the centered discretization. More generally the method applies to systems 

Ut + [u^u,v)]^ = [A{u,v)]^, (10) 
Vt + [v<^{u,v)]^ = [B{u,v)]^, (11) 

which are split into the two subsystems 

ut + [u^u,v)l^0, (12) 
vt + [v^u,v)l^0, (13) 
where $(u, v) plays the role of numerical velocity and 

Ut^[A{u,v)l, (14) 

Vt^[B{u,v)l. (15) 

Systems (12,13) is a family of degenerate systems considered in [2]. In 
particular the scheme in this paper gives neat results for the system (4) in 
reference [2] which models an elastic string problem. 

The numerical scheme. The real line is divided into intervals /j =]ih — 
ih + i E Z. We set t„ = nrh for r small enough. We will construct 
step functions u{x,t), v{x,t) depending on h, which are constant on the 
rectangles /jX]t„, whose step values are denoted respectively. 

The indices h are often skipped to simplify the notation: u stands for Uh, .... 
If a < 6 one sets 

L(a, 6) := length of [0, 1] n [a, b], (16) 

i.e. 

L(a, b) = max{0, min{l, b) — max{0, a)). (17) 

The notation L allows a synthetic formulation of the transport, without being 
forced to distinguish several cases depending on the signs of the numerical ve- 
locities. By induction we assume that the set of values {uf, v^}i^z is known. 
We obtain the set of values {u^~^^,vf~^^}i^z as follows. 

•First step: transport with velocity $ during time rh 

$r :=$«,<), (18) 
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u,:=<_iL(-l+r$ti,r$?_i)+<L(r$^l+r$r)+<+ii^(l+r$r+i,2+r$r^^^ 

(19) 

V, := <_,L(-l+r$t„r<l>ti)+<L(r$r, l+r<l>I')+<+iL(l+r$IVi, 2+r$^+i). 

(20) 

When the CFL condition r|$"| < 1 Vz, Vn is satisfied the first terms in 
(19,20), when multiphed by h, represent the quantities u,v issued from the 
cell li-i between times tn and tn+i that lie in the cell li at time tn+i- Indeed 
the cell /j-i = [(i — |)/i, {i — ^)h] has been transported according to the vector 
r^2_ih, since is the numerical velocity and the duration time is rh. The 
overlap with the fixed cell /j = [{i — |)/i, {i + |)/i] has a length of r$"_]^/i if 
> 0, if < 0, taking into account the CFL condition rl^-'.J < 1. 
Prom (16) one finds L(-l + r<^'^_^, r<l>^_J = r^1_^ if > 0, if < 0. 
Division by h is due to the fact that Ui, Uj are mean values on cells of length h. 

The second terms in (19,20), when multiplied by h, represent the quanti- 
ties u, V issued from the cell /j that remain in /j at time tn+i- Indeed the cell 
[{i — |)/i, {i + |)/i] has been transported by the vector r$"/i. The overlap 
with the fixed cell [{i - |)/i, {i + l)h] is h - r^fh if $f > 0, /i + r^fh if 
$7 < 0. From (16) one finds L(r$f , 1 + ) = 1 - if $f > 0, 1 + 
if < 0. 

The third terms are similar to the first ones: they concern the quantities 
u, V issued from the cell /j+i that lie in the cell /j at time tn+i, with the same 
verification as above. Note that Ui,Vi depend on n, which is not explicitely 
stated to shorten the notation. 

•Averaging step. For a value a, < a < 0.5, to be chosen, we set 

Ui :^ aui-i + {1 — 2a)ui + aui+i, (21) 

Vi -.^ avi^i + {1 — 2a)vi + avi+i. (22) 

In the case A — 0, B — the averaging step is useless. Indeed the idea under- 
lying the elaboration of the scheme is that the first step works well without 
averaging, and that the numerical defects of the centered discretization in 
the last step should be compensated by the averaging step performed before 
it. The splitting should be chosen so as to minimize the importance of the 
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terms involved in the last step. 

•Last step: centered discretization 

:= + - AK-i,vl,)], (23) 

v?^' v^ + ^[^«+i,^lVi) - B{K-i,vU)]. (24) 

The scheme works well for singular shocks and delta shocks. The theorem 
below shows that it gives an approximate solution of the equations. 



3. Statement of the theorem. Let T > be given. Let us seek a 
solution on M X [0,T]. The initial conditions it°,v° are discretized as usual 
by mean values in the cells since they are supposed to be functions. Let 
us apply the scheme under the assumptions (25-29) below: there exists a 
sequence of values /i, /i — ?> 0, a corresponding sequence of values r, r > 0, 
and real numbers /3, 7 G [0, 1[ such that when h ^ 

- ^ (25) 
r 

yn<^Wi r\m < 1, (26) 

~ rh I J I - > V / 

which is the CFL condition, 

yn<^yi /I'^l^^l = 0(1), (27) 
which is a constraint on the numerical velocity allowing it to tend to infinity, 

i i 

which is the L-'^-stabihty in m, v, 
^^^^^^ J]|A«,<)|/i^+^ = 0(l),J]|S«,<)|/i^+^ = 0(l). (29) 



Theorem. Consistence of the scheme. As long as (25-29) are satisfied 
then the scheme is consistent on Mx]0,T[ in the sense of distributions, i.e. 
ifuh,Vh, are the step functions from the scheme, then, V'i/' £ C^{M.x]0,T[), 
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/ + vM, - AMi,M=^dt ^ 0. (30) 

J [vhipt + Vh^{uh, Vh)ipx - B{uh, Vh)ipa]dxdt 0, (31) 
when h ^ 0. More precisely the integrals in (30,31) are equal to 

0{^) + 0{h^-f^) + 0{h^-^). (32) 

The scheme will be of order one in the usual cases in which r is constant, 
/? = 7 = 0, but of an order strictly less than one for singular shocks from the 
fact that the values of the numerical velocity increase when /i ^ 0, which 
forces r ^ and ^ > 0. 

4. Proof of the theorem. / uiptdxdt = J2i,n ^? Leik „ i'tdxdt = J2i,n K[i'4^t)2+ 

0{h)]rh' = E,n <^^^rh' + u20{rh)rh' + <0(/z)r/i^ 
Since | Ei,„ <C>(/i)r/i2| < En^^^Ei K\\0(h)\h < const.T\0(h)\ from (28), 
one obtains 

f ui^tdxdt = - <^^)^C + 0{h). (33) 

i,n 

Similarly 

/ $(«, v)u^P,dxdt = E,. ^^dxdt = E,. $r<(^x)^/iHE,,n ^y70{h)rh' = 

E,n $>r^^^P^r/i2 + E,„ $>70(/.)r/i2 + E,„ $>rO(/i)r/.2. 

IEin^>"C>W^^^I ^ En^^Ei l^rlKI|0(/i)|/i < const.Th-^h < const.h}-^ 
from (27,28). Finally 

j ^{u,v)u^,dxdt = -hY,m< - *r-i<-i)^C + 0{h^-^). (34) 
Similarly 

E,n^«><)0(/^)r/l^ = E,n^«><)%^^/^' + E,n^«.^r)0(/i)r/l^. 
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I j:,,^A{u2,vf)0{h)rh'\ < T.n^hY.i \A{vJl,v2)\\0{h)\h < const.T.h-^h < 
const. h}~'^ from (29). Therefore 

/ A{u, v)i^,dxdt = J2 rh[A{uU, <_i) - <)]C + 0(/i^-^). (35) 

i,n 

Setting 

/ := y [utpt + u^{u, v)'ijjx — A{u, v)ipj^dxdt, (36) 
one finally obtains from (33-36) 

r(A«, <) - + 0{h) + O(/ii-0 + 0(/.i-^'). (37) 



Up to this point the formulas of the scheme have not yet been used. From 
(23,21) 

<+i = + a{u,_, - 2ui-rUi+,) + i[A«+i,i;r+i) - A{uU,vU)\- 
Therefore, from (37) 

/ = 7i + 72 + /3 + 0(/i) + 0{h}-P) + 0(/i^-^) , (38) 

where 

h = -/^ $^[«. - < + r«$,'^ - <_i$r-i)]C, (39) 

i,n 

h = -ha ^(wi-i - 2ui + Mi+i)^f , (40) 

i,n 

^3 = -^EM^(<+i,^m)-^(<-i>^l^i)-2K<,<)-^(<-i,^ti)]}C- 

i,n 

(41) 

We are going to prove successively bounds for Ji, /2, /a. 

• bound for Ii . In Ji fix an index io and consider successively the two 
cases < and > 0. 

If < then, from (16) and the CFL condition (26), L{r^l, l + r$^J = 
1 + L{l + r^l,2 + r^l) = -r^l and L(-l + r*^^, r$^J = . Therefore 
from (19) 
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— u2^{l + r$f^)+terms not involving uf^, 
UiQ-i — — ii^r$"^+terms not involving u^, 
UiQ+i does not involve uf^. 

From the CFL condition the other terms Ui do not involve m"^. Therefore in 
the sum WiV^f the term uf^ occurs in (and only in) 

<o(i+^*"o)C-<o^*roCo-i- 

Consequently in the sum ~ "^i + r{u2^2 — u2_i^f_i)]il)^ , the term 

involving uf^ is 

<o (1 + rK)^l - Kr^li^l-i - <Co + rKoK€o - r^K^o+i (42) 



where the first two terms come from and Ui^-i. The sum (42) is equal 
to ruK^Dlijl - V'fo-i + i^Z - i^Z+i] = r<^$r„0(/i') from Taylor's formula 
applied to ip. 

If > then, an analogous reasoning involving and Ui^+i instead of 
and Uig-i gives the value 0. Therefore from (39) 



\h\ < hZi,,n<r^lO{h') = ZnrhZi^?u2hO{h), i.e. from (27,28) 

l^ = 0{h^~f'). (43) 



• bound for I2. Prom (40) I2 — —ha^-^Ui{iJj'^_^i — 2-?/'^ + ipl 

« En Ei^iO{h^) = «T^0(1) since one has j^i < D K\h = 0(1). 
Indeed (19) implies the formula 

(44) 

The definition (16) of L implies L(— 1 + a, a) + L{a, l + a) + L(l+a, l+2a) = 1. 
Therefore from (44) \ui\ < This implies 

/2 = 0(^). (45) 
. bound for h. h = -f E,,„{^«,0^r-i-^«,<)V'r+i-2A«,<)^,"+ 

from Taylor's formula in ■0 and (29). Therefore 
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Is^O{h'-^). 

Finally from (38,43,45,46) 



(46) 



I = 0{h'-^) + 0{h}-^) + O(^), (47) 
which ends the proof. □ 

5. Approximation of the Keyfitz-Kranzer system (1,2). We con- 
sider successively the three different typical solutions of Riemann problems 
in figures 8,7,6 in [4]: singular shock, intermediate overcompressive shock 
and usual shocks. The numerical solutions obtained from the scheme are 
identical to those shown in [4] even in absence of additional viscosity. We 
consider first the Riemann problem in figure 8 in [4] , which shows a singular 
shock. The initial data is {uuVuUr.Vr) = (1.5,0,-2.065426,1.410639). We 
adopt the values a — 0.2, /3 = 0.5,7 — 0-4- One chooses the value of Vh 
close to the maximum value of r that satisfies the CFL condition (26). For 
simplicity we denote 

"(27)" := h^maXi^nKl, 
"(28)" := maXr.{J2\<\h,J2K\h), 

i i 

"(29)" := maXniY. <)|/.^+^ 01/^^+^) 

i i 

for the values in the assumptions of the theorem. 

In order to check the consistence theorem we present the values of - 
that must tend to from (25), and the values "27", "28", "29" that must be 
bounded. Results of a test for T = 5 with the interval [—4, 4] are given in 
the table below. 



h 


r 


k 


"(27)" 


"(28)" 


"(29)" 


0.0400 


0.300 


0.1333 


0.6289 


14.97 


3.62 


0.0200 


0.240 


0.0833 


0.5830 


14.97 


2.84 


0.0100 


0.170 


0.0588 


0.5309 


14.96 


2.26 


0.0050 


0.132 


0.0379 


0.5271 


14.93 


1.84 


0.0025 


0.095 


0.0263 


0.5178 


14.90 


1.53 


0.00125 


0.065 


0.0192 


0.5021 


14.87 


1.29 


0.00062 


0.040 


0.0156 


0.4326 


14.85 


1.10 


0.00031 


0.025 


0.0125 


0.4024 


14.83 


0.96 
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Now we choose T — 1 and the interval [—0.5, 0.5] in order to reach smaller 
values of h. The values of the parameters are again a — 0.2, /3 — 0.5,^ — OA 



h 


f 


h 

r 


" ('97V" 




"('99V' 


0020 


18 


01 1 1 

\J .\J J. J. J. 


0.2444 


1.9232 


1791 


0.0010 


0.13 


0.0077 


0.2337 


1.9170 


0.1480 


0.0005 


0.09 


0.0056 


0.2225 


1.9109 


0.1252 


0.00025 


0.06 


0.0042 


0.2090 


1.9054 


0.1081 


0.000125 


0.043 


0.0029 


0.2070 


1.8999 


0.0973 


0.0000833 


0.035 


0.0024 


0.2051 


1.8972 


0.0926 


0.0000625 


0.030 


0.0021 


0.2028 


1.8955 


0.0898 


0.0000500 


0.026 


0.0019 


0.1979 


1.8944 


0.0874 


0.0000333 


0.021 


0.0016 


0.1955 


1.8923 


0.0848 


0.0000250 


0.019 


0.0013 


0.2010 


1.8907 


0.0847 


0.0000166 


0.015 


0.0011 


0.1957 


1.8891 


0.0829 


0.0000125 


0.012 


0.0010 


0.1847 


1.8884 


0.0803 



insert figure 1 

Figure 1: The numerical solution from the last test in the second table. 
One can observe that the scheme reproduces exactly the aspect of the singular 
shock in figure 8, in reference [4]. 

The values of r are chosen close to the maximum values for which the 
scheme satisfies the CFL condition r||M||oo < 1 . One observes that the 
quantity - ^ as \/^ and that the three quantities in the columns "27", 
"28" ,"29" are bounded (since quantity "27" is proportional to the sup. of 
\u\ it is very sensitive to the chosen value of r close to the sup. of values 
of r that satisfy the CFL condition). Therefore, since /3 — 0.5,7 — 0-4, the 
scheme is of order 0.5 in h from (32). This is not a good result in general 
from a numerical viewpoint; however the presence of singular shocks gives a 
numerical velocity which is of the order instead of a constant in the usual 
situations in which the scheme is always of order 1. We can see also that 
the bounds in the proof of the theorem are not optimal since one has used 
a bound involving the factor ||$"||oo while $" = "u" is uniformly bounded 
independently of time except on the singular shock. Indeed one can see that 
the scheme gives acceptable results. On a standard PC top values of the 
peak in v in the above tests have reached the value 3700 for the Riemann 
problem under consideration while they have reached values 10^ in the case 
of system (3,4). For the system (3,4) we will rigorously prove in section 6 
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that the scheme is of order 1. The set of results in these two tables gives a 
reasonable presumption that the decrease of ^ and the boundedness of the 
three quantities "27", "28", "29" continue to hold when /i — ?> 0, which would 
allow the theorem to be applied with confidence. If we only consider values 
of h for which (25-29) have been tested then the proof of the theorem gives a 
bound (depending on the sup norm of the derivatives of order two of ip and 
its support) for the integrals in (30,31) according to (32). 

Then we consider the Riemann problem {ui, vi, Ur, Vr) = (1.5, 0, —1.895644, 1.343466) 
represented in figure 7 in [4]. In this case we choose a = 0.2, /3 = 0,7 = 0. 
We obtain the following table: 



h 


r 


h 


"(27)" 


"(28)" 


"(29)" 


0.0050 


0.45 


0.0111 


1.9205 


1.6913 


1.2831 


0.0010 


0.45 


0.0022 


1.9205 


1.6965 


1.2753 


0.0005 


0.45 


0.0011 


1.9205 


1.6972 


1.2743 


0.00025 


0.45 


0.0006 


1.9205 


1.6975 


1.2739 


0.000125 


0.45 


0.0003 


1.9205 


1.6977 


1.2736 


0.0000625 


0.45 


0.0001 


1.9205 


1.6977 


1.2735 



insert figure 2 

Figure 2: The numerical solution from the Riemann problem considered 
in the third table (h — 0.04, r = 0.45^. One can observe that the scheme in 
this paper reproduces exactly the aspect of the limit overcompressive shock in 
figure 7 in [4]. An enlargement has been done in the horizontal direction to 
observe the detailed structure of the shock. 

The results are very clear due to the boundedness of u in this case. There 
is a very natural presumption that these results continue to hold when h ^ 0. 
One can see that the scheme is of order one in h as this follows from the the- 
orem. 

For the third Riemann problem, {ui, vi, Ur, Vr) = (1.5, 0, —1.725862, 1.276293) 
in figure 6 in [4], in which there is no singular shock, the results are very clear, 
exactly the same as those in the above table. We have always observed re- 
sults as good in the case of bounded numerical velocity. 

We now present a system for which a full proof of consistence in the sense 
of distributions has been obtained. 
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6. Application to system (3,4). 

• One considers the 2x2 system (3,4) which produces delta- waves in 
the variable v, see [3]. Here ^{u,v) = u, A = B = 0. In this case one 
can choose a = in (21,22) since the last step (23,24) is absent. Then 
-u"^^ = Ui, v'^'^'^ = Wi, the choice a > works as well with the same proofs. 
It follows from (44) that Xlj — Xlj 1^11- Therefore by induction on 
n 'Ylii ^ Si I'^i'l- The same proof applies for v. Choosing the ini- 

tial condition in this proves (28). To prove (25,26,27) we will prove 

the maximum principle in the numerical velocity u. 

Lemma. If rmaxi\u^\ < | then u satisfies the maximum principle. 

Proof. Let the index i be fixed. Consider the various possible combina- 
tions of signs in the three values u".^, tt", u^j^-^. In each case one will check 
that 

mm«_i,<,<+i) <Ui^ <+i < maa:«_i,<,<+i) 

which proves the maximum principle by induction on n. By induction up to 
order n the condition rmaxi\v!l\ < \ implies rmaxi\u^\ < \. Now we pass to 
order n + 1. 

• case (-|-,-|-,-|-). Formula (19) with ^ — u gives 

ur = uUruU + <(1 - r<) = < + r{uU - <)«-i + <)• (48) 

First note that Tlj > because 1 — ru^ > from the property rmaxi\u^\ < |. 
We consider successively the two cases vl^ > u^_i and vl^ < u^_i. If > "u^.i 
then (48) gives Ui < u^. If < u^_^ then Wj — u'^_^ = (m" — u^_i)\}- — r{u^ + 
M^_i)] < since the last factor is > by induction. We have checked that 

< Mi < max{u1_i,uD. 

• case (-|-,-|-,-). Formula (19) gives 

u, = uUruU + <(1 - r<) + <+i(-r<+i). (49) 

First let us prove that Ui > it^+i. The properties u^_^ > 0, > 0, ruf < \ 
imply that 

Ui > Ui_^_^{—ru'^_^_^) > u'^_^_^ since < —ru^_^_i < | and u'^_^_^ < 0. 
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Now let us check that Ui < max{u^_i^u^). Formula (49) and u^j^^^ < 
imply Ui < u'^_irUi_i + — ru^). Prom this inequality the proof is the 
same as in the case (+++). 

• case (-,+,+). Formula (19) gives Ui = m"(1 — rw") which implies Ui < 
since < ru^ < \ and, rtj > 0. 

• case (-,+,-)■ Formula (19) gives 

Ui = <(1 - r<) + <+i(-r<i) = < + r[-«+i)^ - «)^] < <. 

Now Ui - <+i = < - <+i - r[«+i)2 + (^n)2]_ Sij^^g ^n^r^^ < (^^^^^2 ^ 

«+i)' < «)' + «+i)' - 2«+i = « - <+J^ Therefore Ui - > 

~ ~ ''(■^r ~ "K+if = ~ ■^r+i) [1 ~ ''(■^r ~ ■^r+i)] ^ O since the second 
factor is positive, which implies Ui > wf+i- 

In the four cases in which uf < the verifications are similar. 

Finally we have proved properties (25-28), with r independent oih, /3 — 
and 7 = since A — B = 0. Therefore from the theorem the scheme con- 
verges in the sense of distributions and is of order one in h. It has been 
checked numerically that its real interpretation is the well known solution. 

6. Conclusion. We have presented a numerical scheme which captures 
the singular shock solutions of the Keyfitz-Kranzer model without recourse 
to a vanishing viscosity method. We have observed numerically exactly the 
same results previously observed by the various authors. The consistence of 
the scheme for this system has been checked numerically up to very small val- 
ues of h. The theorem states that the approximate solutions from the scheme 
tend to satisfy the equations in the sense of distributions. This scheme adapts 
to degenerate systems such as the Korchinski model system and the Keyfitz- 
Kranzer system of elasticity. In the case of the Korchinski system consistence 
in the sense of distributions has been fully proved. 
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